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The role of the many-body (spin-polaronic) effects in the scanning tunneling spectroscopy of half- 
metallic ferromagnets (HMF) is considered. It is shown that the non-quasiparticle (NQP) states 
exist in the majority or minority spin gap in the presence of arbitrary external potential and, 
in particular, at the surfaces and interfaces. Energy dependence of the NQP density of states is 
obtained in various models of HMF, an important role of the hybridization nature of the energy gap 
being demonstrated. The corresponding temperature dependence of spin polarization is calculated. 
It is shown that the NQP states result in a sharp bias dependence of the tunneling conductance near 
zero bias. Asymmetry of the NQP states with respect to the Fermi energy provides an opportunity 
to separate phonon and magnon peaks in the inelastic spectroscopy by STM. 



I. INTRODUCTION 



The history of the investigations of half- metallic ferromagnets (HMF) starts from the electronic structure calculation 
for NiMnSb 0]; later a number of other examples were discovered, e.g., CrC>2, Fe3C>4, a number of the Heusler alloys 
Co2MnZ and RMnSb (for a review, see Refs. a a. These substances have metallic electronic structure for one 
spin projection (majority- or minority-spin states), but for the opposite spin direction the Fermi level lies in the 
energy gap. Owing to this fact HMF attract now a growing attention of researchers in connection with "spintronics" , 
or spin-dependent electronics |^|. The spin- up and spin-down contributions to electronic transport properties have 
different orders of magnitude, which can result in a huge magnetoresistance for heterostructures containing HMF 0. 
Some evidences of the HMF behavior in colossal magnetoresistance (CMR) materials like Lai-^Sr^MnOa found by 
using tunneling spectroscopy 0, ||| and photoemission technique have increased considerably the interest in the 
half-metallic ferromagnetism; however, the situation in the CMR systems is controversial, as demonstrate Andreev 
reflection experiments [8j. 

Peculiar band structure of HMF results in an important role of incoherent (non-quasiparticle, NQP) states which 
occur because of correlation effects @. The appearance of NQP states in the energy gap near the Fermi level 

[ToL ITTL I12I is one of the most interesting correlation effects typical for HMF. The origin of these states is connected 
with "spin-polaron" processes: the spin-down low-energy electron excitations, which are forbidden for HMF in the 
one-particle picture, turn out to be possible as superpositions of spin-up electron excitations and virtual magnons. The 
density of these states vanishes at the Fermi level Ep for temperature T — 0, but increases drastically at the energy 
scale of the order of a characteristic magnon frequency ZU. The existence of NQP states is relevant for spin-polarized 
electron spectroscopy 0, , NMR 141, core- level spectra of the HMF 0], and subgap transport in ferromagnet- 
superconductor junctions (Andreev reflection) |l6j . Several experiments could be performed in order to clarify the 
impact of the NQP states on spintronics. In particular, I — V characteristics of half-metallic tunnel junctions for the 
case of antiparallel spins are completely determined by NQP states 0] . Recently the density of NQP states has 
been calculated from first principles for a prototype HMF, NiMnSb |l9j |. and for CrAs pfT ]. 

On the other hand, HMF are very interesting conceptually as a class of materials which may be convenient to treat 
many-body solid state physics that is essentially beyond band theory. It is accepted that usually many-body effects 
lead only to renormalization of the quasiparticle parameters in the sense of Landau's Fermi liquid (FL) theory, the 
electronic liquid being qualitatively similar to the electron gas (see, e.g., Refs l2ltl22j) . On the other hand, NQP states 
in HMF are not described by the FL theory. As an example of highly unusual properties of the NQP states, we note 
that they can contribute to the T-linear term in the electron heat capacity |12ll23|. despite their density at Ep is zero 
for T = 0. 

Spin-polarized scanning tunneling microscopy (STM) p4| is a very efficient new method which enables one to probe 
directly the spectral density with spin resolution in magnetic systems. The spin-polarized STM should be able to 
probe the NQP states via their contribution to the differential tunneling conductivity dX/dV . At zero temperature, 
NQP states arise only above Ep for the case of minority-spin gap and only below Ep for the majority-gap HMF 
0. Unlike the photoemission spectroscopy which probes only occupied electron states, STM detects the states both 
above and below Ep, depending on the sign of bias. 

Theoretical investigation of NQP contributions to STM spectra is the aim of the present paper. The paper is 
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organized as follows. In Sect. 2 we discuss a general expression for the tunneling current X as applied to HMF. In 
Sect. 3 the effect of surface potential and other spatial inhomogcncitics on the NQP spectral density is considered. In 
Sects. 4 and 5 we calculate the energy and temperature dependences of dl/dV and treat the problem of tunneling- 
current spin-polarization at finite temperatures. In Sect. 6 the bias dependences of the tunneling conductance are 
discussed. 



II. CALCULATION OF THE TUNNELING CURRENT 



A general expression for the tunneling current in the lowest order in the tunneling matrix elements has the form 

mug 
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where e is electron charge, a is the spin projection, V is the bias, / (e) is the Fermi distribution function, Greek 
(Latin) indices label electron eigenstates for the sample (tip) Vw,?/w 



M ™ = £ / dA «* V ^ ~ CVi„) (2) 

is the current matrix element, m is the free electron mass (the surface integral in Eq. is taken over arbitary area 
between the tip and the sample), and 

Al (E) = -hmG* (E) (3) 
is the corresponding spectral density, G° (E) = G° v (E) is the sample Green's function, 

G? vX (E) = {{c va \c{ a )) El 

cJ, CT being the creation operators for conduction electrons. It is worthwhile to emphasize that the expression JQ) takes 
into account correlation effects in a formally exact way, assuming that the tunneling probability is small. In fact, 
the latter condition should be satisfied for proper STM measurements, otherwise they cannot be considered as a true 
probe. In the WKB approximation Eq.JTJ takes the form j2]j 
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where is the average of sample and tip work functions (which is assumed to be large in comparison with eV for 
simplicity), z is the distance between the surface and the tip. Here N t (E) is the density of states (DOS) of the tip 
material, which is supposed to be smooth and thus is replaced by its value at the Fermi energy, and 

g°{E) =Y,K (*0M * (k||)K> (5) 

V 

is the density of states of the sample with zero in-plane component of the wave- vector: k|| = so that the summation 
is performed only over two points of the Fermi surface. This condition means that the tunneling probability has 
a sharp (at not too small z) maximum for the states with velocity direction normal to the surface. For a generic 
multi-sheet Fermi surfaces the condition of the dominant tunneling is, generally speaking, more complicated, but this 
modifies only some weakly bias-dependent factors. 



III. NON-QUASIPARTICLE STATES IN INHOMOGENEOUS MATERIALS 

Since STM probes only surface one has to discuss first modification of NQP states in comparison with the case of 
ideal bulk crystal. The existence of NQP states at the surface of HMF has been demonstrated for a narrow-band 
Hubbard model |2S| . Here we present a general derivation valid in the case of arbitrary inhomogeneity. 
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To describe the effects of electron-magnon interaction for the inhomogeneous case we use the formalism of the 
exact eigenfunctions developed earlier for the impurity-state problem in a ferromagnetic semiconductor 29] . The 
corresponding Hamiltonian of the s — d exchange model reads 

K = -^V 2 + [/ CT (r) (6) 
2m 

where U a (r) is the potential energy (with account of the electron-electron interaction in the mean field approximation) 
which is supposed to be spin dependent, x P cr (r) is the field operator for the spin projection cr, S(r) is the spin density of 
the localized-moment system, <5S(r) = S(r)— (S(r)) is its fluctuating part, the effect of the average spin polarization 
(S(r)) being included into U a (r). We use an approximation of contact electron-magnon interaction described by the 
s — d exchange parameter /, 

TC d = - J q S q S_ q (7) 
q 

is the Heisenberg Hamiltonian of localized spin (for simplicity we neglect the inhomogeneity effects for the magnon 
subsystem) . 

Passing to the representation of the exact eigenfunctions of the Hamiltonian Hq , 
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P 

* CT (r) = ^Vw(r)<W, (8) 



one can rewrite the Hamiltonian 10) in the following form: 

Ti = ^ e vcrc\ a c va -I (i/a, fJ,/3\q) 5Sqcl a cr a 0C^ + H d (9) 

where 

{va,na' |q) = {fia'\ e lqr \va) . 

We take into account again the electron-spectrum spin splitting in the mean-field approximation by keeping the 
dependence of the eigenfunctions on the spin projection. 

We restrict ourselves to the spin- wave region where we can use for the spin operators the magnon (e.g., Dyson- 
Maleev) representation. Then we have for the one-electron Green's function 

Gl(E) = [E-e„ t7 -K{E)r\ (10) 

with the self-energy S^(-B) describing correlation effects. 

We start with the perturbation expansion in the electron-magnon interaction. To second order in / one has 

K(E) = 2I 2 SQ:(E) (11) 

with 

QUE) =X)IM,/U |q) l 2 "**"} ME) = £ | {v 4, (J, T |q) f 1 + N «-< (12) 

where = /(e^cr) , w q is the magnon energy, iV q = JVfl(w q ) is the Bose function. 
Using the expansion of the Dyson equation ijTUI) we obtain for the spectral density 

A V *(E) = --lmGZ(E) =5(E-E va ) 

-6\E-e va )RzEZ{E)-- » K > (13) 

7T yH/ &ua) 
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The second term in the right-hand side of Ea. (|13fl gives the shift of quasiparticle energies. The third term, which arises 
from the branch cut of the self-energy, describes the incoherent (non-quasiparticle) contribution owing to scattering 
by magnons. One can see that this does not vanish in the energy region, corresponding to the "alien" spin subband 
with the opposite projection — a. 

Neglecting temporarily in Ea. H12|) the magnon energy cj q in comparison with typical electron energies and using 
the identities 

y\iyMt n = _1 f^FW y. |Hq)P 

^ E-e u v ^' 7T J E-E' ^E'-s u + iO v ' 

J q 

= -- I dE'P^l^{ v \{E' -Ho+iO)- 1 ^) 
7r j tL — tj 

we derive 

Ej(£) - 2I 2 S J dE'f(E') (u ]\5(E-E'- H l ) \v \) (15) 
Zl(E) = 2I 2 S JdE[l- f(E')} {v l\s(E-E'- H T ) \v j) (16) 

Here we restrict ourselves only to the case of zero temperature T = neglecting the magnon occupation numbers -/V q . 
Using the tight-binding model for the ideal-crystal Hamiltonian we find in the real-space representation 

5^,(25) = 2I 2 S J dE'f(E') f-ilmGj^OE')) <5 R , R , (17) 

X RR ,(E) = 2I 2 S J dE' [1 - f(E')] f-hmGl^E')) 6 R<n , (18) 

where R, R' are lattice site indices, and therefore 

K(E) = J2\^( K )\ 2 ^kn(E). (19) 
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Following the method developed by us earlier 12, 29] one can generalize the above results to the case of arbitrary 
s — d exchange parameter. Simplifying the sequence of equations of motion (cf. Ref l29f) we obtain the integral equation 

(E - e v „)GZ x (E) = 6 vX + aIRZ x (E) - al^E - e K - a )K x (E) GJ„ (E) (20) 

K 

where 

Rlx(E) = E^i' A tl-q)(^t,Ml|q) 
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E - Efj.i + UJq 

Rtx(E) = Yl 0* T, A i |q) {v |, M T |q) „ ^— ^ (21) 

Note that the equation 1)20(1 is exact in the case of empty conduction band (n vcr = 0, one current carrier, ferromag- 
netic semiconductor situation), and for finite band filling this corresponds to a ladder approximation in the diagram 
approach. 

Similar to Q14[l. we obtain after neglecting the magnon energies in (|21|l the equation for the Green' function 

J2(E - e K ^)R: x (E) GZ k (E) = (H R° (E) (E - HI + iD)-^ (E) \ua) (22) 



where we use the matrix notations. Then we have for the operator Green' function 



G a (E) = 



E-Hz + „mz-H- Q n 1+a l Rtr{E) R°{E) 



(23) 
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If we consider spin dependence of electron spectrum in the simplest rigid-splitting approximation e U(J — e v — a I (5 Z ) 
and thus neglect spin-dependence of the eigenfunctions Vw (R-) the expressions (I15H . (|16I) are drastically simplified. 
Then the self-energy does not depend on v and we have 



T, a (E) = 



2I 2 SR a (E) 
l + aIR a (EY 



R\E) 



E - e mT 



(24) 
(25) 



If TLq is just the crystal Hamiltonian (y — k, e va — tku, iko- being the band energy), the expression (12,'SH coincides 
with that obtained in Ref. ^3 f° r the Hubbard model after the replacement J — ► U. 
The expression (|23[) can be also represented in the form 



G a (E) = 



E - Ho* - {K - Ho 1 



l + aIR°(E) 



(26) 



The equation Q26[) is convenient in the narrow-band case. In this limit where spin splitting is large in comparison with 
the bandwidth of conduction electrons we have nl-Hl = -21 5 and we obtain for the "lower" spin subband with 
(7 = —sign/ 

G a (E) = 

For a periodic crystal Eq. (|27|l takes the form 



e 
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R°{E) 



(27) 



This expression yields exact result in the limit / 



E — tk- 



-oo, 



Gi (E) 
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R(e) 



25 
R a {E) 



(28) 



(29) 



with e = E + IS, the bare electron spectrum. In the limit / — > — oo Ea, 128|l gives correctly the spectrum of 
spin-down quasiparticles, 



Gi (E) = 



25 



25+1 



\e-tt 



(30) 



with e = E — I(S + 1), t£ = [25/ (25 + l)]tk- However, it does not describe the NQP states quite correctly, so that 
more accurate expressions can be obtained by using the atomic representation |30| . 



GiiE) 



25 



25+1 



25 
R*(e) 



,R*(e) = 



4^ e - tt 



(31) 



On the other hand, the result 127fl . (|28|l yields a good interpolation description in the Hubbard model [9llTolll2]|. 

The Green's functions (|28|) . 129|) . 131|l have no poles, at least for small current carrier concentration, and the whole 
spectral weight of minority states is provided by the branch cut (non-quasiparticle states) 0, . For surface states 
this result was obtained in Refl28l in a narrow-band Hubbard model. Now we see that this result can be derived in an 
arbitrary inhomogeneous case. For a HMF with the gap in the minority spin subband NQP states occur above the 
Fermi level, and for the gap in the majority spin subband below the Fermi level. 

In the absence of spin dynamics (i.e., neglecting the magnon frequencies) the NQP density of states has a jump at 
the Fermi level. However, the magnon frequencies can be restored in the final result, in analogy with the case of ideal 
crystal, which leads to a smearing of the jump on the energy scale of a characteristic magnon energy ZU. It should 
be mentioned once more that we restrict ourselves to the case of the usual three-dimensional magnon spectrum and 
do not consider the influence of surface states on the spin-wave subsystem. The expressions obtained enable us to 
investigate the energy dependence of the spectral density. 
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IV. THE NON-QUASIPARTICLE DENSITY OF STATES 



An analysis of the electron-spin coupling yields different pictures for two possible signs of the s — d exchange 
parameter /. For I > the spin-down NQP scattering states form a "tail" of the upper spin-down band, which starts 
from Ep (Fig.l) since the Pauli principle prevents electron scattering into occupied states. 

For I < spin-up NQP states are present below the Fermi level as an isolated region (Fig. 2): occupied states with 
the total spin 5 — 1 are a superposition of the states |5}| j) and \S — 1)| T). The entanglement of the states of electron 
and spin subsystems which is necessary to form the NQP states is a purely quantum effect formally disappearing 
at S — > oo. To understand better why the NQP states are formed only below Ep in this case we can treat the 
limit I = — oo. T hen the current carrier is really a many-body state of the occupied site as a whole with total spin 
5 — 1/2, which propagates in the ferromagnetic medium with spin 5 at any site. The fractions of the states |5)| |) 
and |5 — 1}| t) in the current carrier state are 1/(25 + 1) and 2S/(2S + 1), respectively, so that the first number 
is just a spectral weight of occupied spin-up electron NQP states. At the same time, the density of empty states is 
measured by the number of electrons with a given spin projection which one can add to the system. It is obvious 
that one cannot put any spin- up electrons in the spin- up site at I = — oo. Therefore the density of NQP states should 
vanish above Ep. 

It is worthwhile to note that in the most of known HMF the gap exists for minority-spin states . This is similar 
to the case I > 0, so that the NQP states should arise above the Fermi energy. For exceptional cases with the 
majority-spin gap such as a double perovskite Sr2FeMo06 [HI one should expect formation of the NQP states below 
the Fermi energy. 

It has been proven in the previous section that the presence of space inhomogeneity (e.g., surface) does not change 
qualitatively the spectral density picture, except smooth matrix elements. Therefore further in this section we will 
consider, for simplicity, the case of clean infinite crystal; all the temperature and energy dependences of the spectral 
density will be basically the same for the surface layer. 

The second term in the right-hand side of Eq. (|13[) describes the renormalization of quasiparticle energies. The 
third term, which arises from the branch cut of the self-energy Y, ua {E), describes the incoherent (non-quasiparticle) 
contribution owing to scattering by magnons. One can see that this does not vanish in the energy region, corresponding 
to the "alien" spin subband with the opposite projection — a. Further on we perform for definiteness concrete 
calculations in the case I > (the case I < differs, roughly speaking, by a particle-hole transformation). Summing 
up Eq. @afl| to find the total DOS N a (E) and neglecting the quasiparticle shift we obtain 



N^E) = J2 



1 - 



2/ 2 57V q 



kq L (*k+qi - *kt) 2 



<*CE-tk T ) 



Nl (E) = 2/ 2 5V 1 + N « " kT 5(E t kT - «,) (32) 

^ (ik+q-L - ^kT - 

where we consider for simplicity only second-order perturbation expression. Basing on a general consideration in the 
previous section one can prove that, actually, this expression holds for arbitrary 7, at least, in the framework of 1/25 
expansion. 

The T 3 / 2 -dependence of the magnon contribution to the residue of the Green's function, i.e. of the effective electron 
mass in the lower spin subband, and an increase with temperature of the incoherent tail from the upper spin subband 
result in a strong temperature dependence of partial densities of states N a (E), the corrections being of opposite sign. 
At the same time, the temperature shift of the band edge for the quasiparticle states is proportional to T 5 / 2 rather 
than to magnetization [ToL |29| . 

The exact solution in the atomic limit (for one conduction electron), which is valid not only in spin- wave region, 
but for arbitrary temperatures, reads ^sj^ 

S+l + ajS*) 1 S-*{S*) 1 

° [ >~ 25 + 1 E + IS + 25+1 E-I(S + lY {M) 

In this case the energy levels are not temperature dependent at all, whereas the residues are strongly temperature 
dependent via the magnetization. 

Now we consider the case T = for a finite band filling. The picture of N(E) in HMF (or degenerate ferromagnetic 
semiconductor) demonstrates strong energy dependence near the Fermi level (Figs. 1,2). If we neglect magnon 
frequencies in the denominators of Ea. H32|) . the partial density of incoherent states should occur by a jump above or 
below the Fermi energy Ep for I > and / < respectively owing to the Fermi distribution functions. An account 
of finite magnon frequencies a-> q = T>q 2 (T> is the spin wave stiffness constant) leads to smearing of these singularities, 
N(Ep) being equal to zero. For \E — Ep\ «wwe obtain 
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N- a {E) 
N a (E) 



1 

2S 



E — Ep 



3/2 



(a(E — Ep)), a — sgnJ 



(34) 



(a = ± corresponds to the spin projections f, J,). With increasing E 1 — Ep\, N- a /N a tends to a constant value which 
is of order of I 2 within the perturbation theory. 

In the strong coupling limit where |/| — * oo we have from (|32[1 



N- a (E) 
N a {E) 



2S 



(a(E - E F )),\E - E F \ >U 



(35) 



In fact, this expression is valid only in the framework of the 1/25-expansion, and in the narrow-band quantum case we 
have to use more exact expressions (|29|I . (I31|) . In numerical calculations, we follow to Ref|30j and smear the resolvents, 



R{E) -► R{E) = J duK{uj)R(E±uj) 



We use the semielliptic magnon DOS K(uj) which is proportional (with the corresponding shift) to the bare electron 
DOS, the maximum magnon frequency being determined by the electron concentration c |3fij . This approximation 
provides the correct behavior near the Fermi level (cf. Ref|33J) , although gives an unphysical shift of the band bottom 
by the maximum magnon frequency. 

The results of calculations are shown in Figs. 3, 4. One can see that in the model with / — > — oo (for 5 = 1/2 this 
is equivalent to the Hubbard model with the replacement <k — * see Refs lTol IT2T1 the "Kondo" peaks |3^| modify 
considerably the picture. Note that the function — (l/7r)ImR*(E), which yields DOS in the lowest-order approximation 
in the electron concentration, does not have such peaks. 

To investigate details of the energy dependence of N (E) in the broad-band case we assume the simplest isotropic 
approximation for the majority-spin electrons, 



i k j — Ep = £k = 



kp 



2m* 



(36) 



Provided that we use the rigid splitting approximation t k j — + A (A = 21 S, I > 0), the half- metallic situation (or, 
more precisely, the situation of degenerate ferromagnetic semiconductor) takes place for A > Ep. Then qualitatively 
the equation H34[) works to accuracy of a prefactor. It is worthwhile to note that, formally speaking, the NQP 
contribution to DOS occurs also for an "usual" metal where A < Ep. In the case of small A there is a crossover 
energy (or temperature) scale 



T* =V(m*A/k F y 

which is the magnon energy at the boundary of Stoner continuum, T* 
equation i|32|) for the NQP contribution reads 



(37) 

lJ{A/E f ) 2 < uj. At \E - E F \ <C 57 the 



SN^E) a 



In 



l + y/{E- E F ) IT* 



\-^{E-Ep)/T* 



-^{E-E F )/T* 



6{E - E F ). 



(38) 



At \E — E F \ -c T* this gives the same results as above. However, at T* -C \E — E F \ <C w this contribution is 
proportional to —^J(E — Ep) /T* and is negative (of course, the total DOS is always positive). This demonstrates 
that one should be very careful when discussing the NQP states for the systems which are not half- metallic. 

The model of rigid spin splitting used above is in fact not applicable for the real HMF where the gap has a 
hybridization origin 0,0. The simplest model for HMF is as follows: a "normal" metallic spectrum for majority 
electrons Ij36|l and the hybridization gap for minority ones, 



ikj. - E f 
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■sgn(£ k )\/£ k 



(39) 



Here we assume for simplicity that the Fermi energy lies exactly in the middle of the hybridization gap (otherwise one 
needs to shift £k — > £ k + Eq — Ep in the last equation, Eg being the middle of the gap). One can replace in Ea. H32|) 
£k+q by v k q, v k = k/m*. First, we integrate over the angle between the vectors k and q. It is easy to calculate 



1 



vpqA 



X 3 - (X 2 



\3/2 



+ x 



(40) 
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where angular brackets stand for the average over the angles of the vector k, X = kpq/m* A. Here we do have the 
crossover with the energy scale T* which can be small for small enough hybridization gap. For example, in NiMnSb 
the conduction band width is about 5 eV and the distance from the Fermi level to the nearest gap edge (i.e. indirect 
energy gap which is proportional to A 2 ) is smaller than 0.5 eV, so that (A/Ep) 2 < 0.1. 
For the case < E — Ep <^oj one has 

The function b(x) is shown in Fig. 5. Thus the behavior N±(E) oc (E — Ep) 3 ^ 2 takes place only for very small excitation 
energies E — Ep <C T* , whereas in a broad interval T* <C E — Ep <C To one has the linear dependence N±(E) oc E — Ep. 



,,5/2 



(1 + 1/) 



5/2 



+ y + y 



3/2 



•r. 



1/ » 1 



(41) 



V. THE TEMPERATURE DEPENDENCE OF SPIN POLARIZATION 

Simple qualitative considerations |34|, as well as direct Green's functions calculations [ill l35j| for magnetic semi- 
conductors, demonstrate that spin polarization of conduction electrons in the spin- wave region is proportional to 
magnetization 

Nt - N\ 
iV T + iV| 

A weak ground-state depolarization 1 — Pq occurs in the case where I > 0. The behavior P(T) ~ (S z ) is qualitatively 
confirmed by experimental data on field emission from ferromagnetic semiconductors 36] and transport properties of 
half-metallic Heusler alloys |37j . 

An attempt was used 38] to generalize the result l|42() on the HMF case (in fact, using qualitative arguments which 
are valid only in the atomic limit, see Eq.|j23J)). However, we will demonstrate that the situation for HMF is more 
complicated. 

In this section we focus on the magnon contribution to DOS l|32|) and calculate the function 

A = E n V^ Ef ^ ^ 

Using the parabolic electron spectrum = k 2 /2m* and averaging over the angles of the vector k we obtain 

. 2I 2 Sm 2 ^ iV q 



^ (?*) - 9 



where p — N^(Ep,T = 0), we have used the condition q <C kp, q* = m*A/kp — A/vp, where A = 2 |/| S is the spin 
splitting. In the ferromagnetic semiconductor we have, in agreement with the qualitative considerations presented 
above: 



25 ' \T C J 

Further on we consider the spectrum model (|36|l . (|39|l where the gap has a hybridization origin. At T <C T* we 
reproduce the result (|45|l which is actually universal for this temperature region. At T* <C T <C lo we derive 

A = S 2/^( & )^ K9 - S ^ K ^Tl4 (46) 

This result distinguishes HMF like the Heusler alloys from ferromagnetic semiconductors and narrow-band saturated 
ferromagnets. In the narrow-band case the spin polarization follows the magnetization up to the Curie temperature 

Tc ___ 

For finite temperatures the density of NQP states at the Fermi energy is proportional to [ll], l23l l34| 

N(E P ) K r dw . X}") (47) 
J smh(w/T) 
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Generally, for temperatures which are comparable with the Curie temperature Tq there are no essential difference 
between half-metallic and "ordinary" ferromagnets since the gap is filled. The corresponding symmetry analysis is 
performed in Ref. for a model of conduction electrons interacting with "pseudospin" excitations in ferroelectric 
semiconductors. The symmetrical (with respect to Ep) part of N(E) in the gap can be attributed to smearing of 
electron states by electron-magnon scattering; the asymmetrical ( "Kondo-like" ) one is the density of NQP states owing 
to the Fermi distribution function. Note that this filling of the gap is very important for possible applications of HMF 
in spintronics: they really have some advantages only provided that T <C Tc- Since a single-particle Stoner-like theory 
leads to much less restrictive (but unfortunately completely wrong) inequality T <C A, the many-body treatment of 
the spin-polarization problem (inclusion of collective spin-wave excitations) is required. 



Now we consider an application of the results obtained above to the tunneling spectroscopy problem. The formulas 
of Sect. 4 for the energy dependence of NQP contributions are, strictly speaking, derived for the usual one-electron 
density of states at Ep, which is observed, say, in photoemission measurements. However, the factor of g° (E), which 
is present in the expression for the tunneling current (J2J , does not influence the temperature dependence, and therefore 
these results are valid for spin polarization from tunneling conductance at zero bias in STM. 

The only difference in the NQP contributions to g° (E) and N a (E) is in that after summation over the magnon 
wavevector q the integration is performed over not in the the whole Fermi surface, but its two points (see Eq.©). For 
a spherical Fermi surface for majority electrons the results differ by the constant factor of the Fermi surface diameter. 
However, the energy and temperature dependences should be the same in a more general case. 

Consider the bias dependence of the tunneling current for zero temperature. One can see from Eq.Q that 



Again, the last proportionality can be strictly justified in the case of a spherical Fermi surface only, but is qualitatively 
valid for arbitrary electron spectrum. 

One should keep in mind that sometimes the surface of HMF is not half-metallic; in particular, this is the case of a 
prototype HMF, NiMnSb [39j. In such a situation, the tunneling current for minority electrons is due to the surface 
states only. However, the NQP states can be still visible in the tunneling current via the hybridization of the bulk 
states with the surface one. The hybridization lead to the Fano antiresonance picture which is usually observed in 
STM investigations of the Kondo effect at metallic surfaces (see, e.g., Refs. l4Cl l4lt 142) . In these cases the tunneling 
conductance will be proportional to a mixture of N a (eV) and L a (eV), L a (E) being the real part of the on-site 
Green's function, 



(V stands for principal value, E is referred to the Fermi energy). Surprisingly, in this case the effect of NQP states 
on the tunneling current can be even more pronounced in comparison with the ideal crystal. The reason is that the 
analytical continuation of the jump in N a (E) is logarithm; both singularities are cut at the energy ZU; nevertheless, 
the energy dependence of L a (E) can be pronounced, see Fig. 6. This is similar to the effect of enhancement of the 
NQP contribution to the x-ray absorption and emission spectra, which was predicted in Ref. Il5l 

Now we discuss in more detail the energy dependence for \E\ <^U. The analytical continuation of the E 3 ^ 2 9(E)- 
contribution to N a (E) yields the contribution (— E) z / 2 9{— E) in L a (E) which is non-zero on the other side with 
respect to Ep (a situation that is formally similar to the electronic topological transition, see Ref 1461) . The one-sided 
linear dependence in N a (E) according to Ea. (|41() corresponds to £Tn \E\ in L a (E). 

STM measurements of electron DOS give also an opportunity to probe bosonic excitations interacting with the 
conduction electrons. Due to electron-phonon coupling, the derivative dN a (E) /dE and thus d 2 I a (V) /dV 2 at eV = E 
have peaks at the energies E — ±oji corresponding to the peaks in the phonon DOS. According to our results (see, 
e.g., Eg. (|32fl . the same effect should be observable for the case of electron-magnon interaction. However, in the latter 
case these peaks are essentially asymmetric with respect to the Fermi energy (zero bias) due to asymmetry of the 
non-quasiparticle contributions. This asymmetry can be used to distinguish phonon and magnon peaks. 



VI. BIAS DEPENDENCE OF THE TUNNELING CONDUCTANCE 



dZ a {V) 
dV 



oc gl (eV) oc N a (eV) 



(48) 




(49) 
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VII. CONCLUSIONS 

In the present paper we have demonstrated that non-quasiparticle states in half-metallic ferromangnets exist not 
only for an ideal crystal, but also in the presence of an arbitrary external potential. In particular, they occur at the 
surface of the half-metallic ferromagnets. These states can be probed by the STM both directly and via their effect on 
the surface states (the Fano antiresonance case). Therefore, they can be observable even for the situation of surface 
"dead layers" where the surface is not half- metallic. The expressions obtained can be used for realistic electronic 
structure calculations of NQP contributions to the electron energy spectrum of the surfaces of HMF. 

Temperature dependence of the spin polarization at the Fermi energy which can be also probed by the STM follows 
the temperature dependence of magnetization an very low temperatures. For the HMF with a hybridization gap, 
there is a crossover energy (temperature) T* <C Tc where the character of the temperature dependence is changed. 
The energy dependence of the NQP contributions (and consequently the bias dependence of the tunneling current) 
is strongly influenced by the band structure too. In particular, for HMF with a hybridization gap this demonstrates 
a linear rather than E 3 ^ 2 behavior in a wide interval. In the narrow-band case a Kondo-like peak (Fig. 4) near the 
Fermi level should be observed in tunneling experiments. 

Due to asymmetry of NQP states with respect to the Fermi energy, the magnon peaks in d 2 T a (V) /dV 2 are also 
asymmetric with respect to the zero bias, in contrast with the phonon ones. This gives an opportunity to distinguish 
between phonon and magnon peaks in the inelastic spectroscopy by STM. 

In principle, the NQP effects discussed should exist also in usual metallic ferromagnets. However, only in HMF 
they can be picked up in a pure form. 

The research described was supported in part by Grant No. 747.2003.2 (Support of Scientific Schools) from the 
Russian Basic Research Foundation and by the Netherlands Organization for Scientific Research (Grant NWO 
047.016.005). 
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FIG. 1: Density of states in the s — d exchange model of a half-metallic ferromagnet with S = 1/2, 1 — 0.3 for the semielliptic 
bare band with the width of W = 2. The Fermi energy calculated from the band bottom is 0.15 (the energy is referred to Ef). 
The magnon band is also assumed semielliptic with the width of u> max = 0.02. The non-quasiparticle tail of the spin-down 
subband (lower half of the figure) occurs above the Fermi level. The corresponding picture for the empty conduction band is 
shown by dashed line; the short-dashed line corresponds to the mean-field approximation. 



N + (E) 




FIG. 2: Density of states in a half-metallic ferromagnet with / = —0.3 < 0, other parameters being the same as in Fig.l. The 
spin-down subband (lower half of the figure) nearly coincides with the bare band shifted by IS. Non-quasiparticle states in the 
spin-up subbands (upper half of the figure) occur below the Fermi level; the short-dashed line corresponds to the mean-field 
approximation. 




FIG. 3: Density of states in a half-metallic ferromagnet in the s-d model with I — ► +oo, S = 1/2. The Fermi energy calculated 
from the bare band bottom is 0.1 (concentration of conduction electrons is c = 0.019). The dashed line shows the function 
-(l/7r)ImR(E). 
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FIG. 4: Density of states in a half-metallic ferromagnet in the s-d model with I — ► —00, S = 1/2. The Fermi energy calculated 
from the bare band bottom is 0.2 (c = 0.034). The dashed line shows the function — (l/7r)ImR*(E). 
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FIG. 5: Plot of the function b(x). 
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FIG. 6: Plot of the imaginary (upper line) and real (lower line) parts of the Green's function near the Fermi level in a 
half- metallic ferromagnet with the same parameters as in Fig.l. 



